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We study the one-dimensional spin-1/2 model with nearest and next-to- nearest-neighbor couplings 
exposed to a homogeneous magnetic field /13 and a dimer field with period q and strength 5. The 
latter generates a magnetization plateau at M = (1 — q/n)/2, which evolves with strength S of 
the perturbation as <5 e , where e = e(/i3,a) is related to the 77-exponent which describes the critical 
behavior of the dimer structure factor, if the perturbation is switched of (5 = 0). We also discuss 
the appearance of magnetization plateaus in ladder systems with / legs. 



(N 

i 

CO 

c3 



1 

C 

O 
o 



> 

(N 
(N 

o 

00 
c3 



c 

o 
o 



X 



I. INTRODUCTION 

The existence or nonexistence of a gap between the 
energies of the ground state and the low lying excited 
states is the most important criterium for the cjiticality 
of a quantum spin system. Haldane's conjectureoa states 
that one-dimensional (ID) quantum spin systems have no 
gap for half integer spin s, but do have a gap for integer 
spin s. The conjecture only holds for an appropriate 
choice of the couplings of the spins at nearest neighbor 
sites.! rjn 

The Lieb, Schultz, Mattis (LSM) constructionfi'El al- 
lows rigorous statements on the degeneracy of the ground 
state. Starting from the unitary operator 



U 



exp 




(1.1) 



it is straightforward to prove that the application of the 
operator U fe , k = 1,2, ... on the ground state |0) gener- 
ates new states 

|fc)=U fc |0) k = 1, 2, k finite, (1.2) 

with an energy expectation value 

(fr|H|fc) - (0|H|0) = 0(JV _1 ) (1.3) 

approaching the ground state energy Eq = (0|H|0) in the 
thermodynamical limit N — > 00. 

Of course the crucial question is whether the new states 
|fc) arc different from the ground state |0) or not. This 
question can be answered by an analysis of the quantum 
numbers of the states \k), k = 1, 2, . . . . For example in 
the case of the Spin-1/2 Hamiltonian with nearest neigh- 
bor couplings: 



H(A 3 ) 



and 



N 

1=1 



N 



2/i 3 S 3 (0), 



S (g) = Xy»'SF, a =1,2, 3, 



(1.4) 



(1.5) 



the ground state has momentum p s — 0, tt and the total 
Spin Sj, = S = NM(h,3), where M is the magnetiza- 
tion. The new states \k) turn out to be eigenstates of the 
translation operator T: 



where! 



T|fc) = TU fc |0> = e lpk \k), 



Pk =Ps + kq 3 (M), 



(1.6) 



(1.7) 



and g 3 (M) = 7r(l - 2M). For M = the ground state |0) 
and the new state |1) differ in their momenta by 7r and 
are therefore orthogonal to each other. 
For M = 1/4 one finds a fourfold degeneracy of the 
ground state with momenta pfe = p s + kn/2, k = 0,1,2,3. 

It should be noted that the LSM construction allows 
to identify the zero frequency e xcit ations (soft modes) in 
the model with Hamiltonian fll.4p . Some of these soft 
modes induce characteristic signatures, e.g. zeroes in the 
dispersion curve and singularities in the transverse and 
longitudinal structure factors at the soft mode momenta 
q = q\(M) = 7r, q = (73 (M), which, can be easily recog- 
nized even on rather small systems.13 Following conformal 
field theory the corresponding critical 77-exponents can be 
determined from the finite-size behavior of the dispersion 
curve at the soft mode momentaJj It is known that the 
^-exponents of the soft mode - i.e. the M dependence of 
773 (M), rji(M) - changes Jj if we add further couplings to 
the Hamiltonian ( jl.4| ). In some cases (see the discussion 
below) the soft mode might disappear completely and a 
gap opens between the states, which were gapless before 
switching on the perturbation. 

The following cases (A.-E.) have been studied so far. 



A. A transverse staggered field 



A gap was foundllHil3 in a transverse staggered field of 
strength hiSi(n), 



1=1 



n{h 3 ,hi) = n(h 3 ) + 2hiSi(n), 



(1.8) 



1 



between the states which differ in their momenta by tt. 
Indeed the operator Si(7r) is invariant only under trans- 
lations T 2 and the eigenstates of the Hamiltonian are 
only eigenstates of T 2 . 

In the free field case (h 3 = 0) the T 2 quantum numbers 
of the ground state |0) and of the LSM state |1) = U|0) 
are the same and the twofold degeneracy of the ground 
state is lifted by the explicit breaking of translation in- 
variance. 

The fourfold degeneracy with momenta p = 0, tt, zhrr/2 
which occurs at M = 1/4 and h\ = is lifted in the 
following manner. The states with p — and p — n are 
even with respect to T 2 . The same holds for the ground 
state |0), which is a linear combination of p — and p = tt 
components. A gap opens to the second state, which is 
even under T 2 . The gap evolves with the strength hi of 
the perturbation as h\. The exponent e = e±(h 3 ) is given 
by the exponent 771 (M) 



added to Hamiltonian (1.4): 



e 1 (h 3 )=2[4- m (M(h 3 ))Y 



(1.9) 



associated with the divergence of the transverse structure 
factor at q = tt. 



The LSM construction with the operator (1.1) leads to 
a second state |1) = U|0) which is degenerate with the 
ground state |0) and which is odd under T 2 . This state 
can be constructed as a linear combination of momentum 
eigenstates with p — ±tt/2. 



B. A longitudinal periodic field 

A longitudinal periodic field S 3 (q) of strength 2h q 

H(h 3 ,h q ) = U(h 3 ) + 2h q S 3 (q), (1.10) 

with 

S 3 (g) = [S 3 (?) + S 3 (-«)]/2, (1.11) 

induces a plateau in the magnetization curve M — M{h 3 ) 
at M = (1 — q/-K)/2, i.e. q has to meet the soft mode 
momentum q = q 3 {M). The difference of the upper and 
lower critical field: 



(1.12) 



A(h q ,h 3 ) = h^-h l 3 ^hf hs \ 



evolves with an exponent, which is again related via ( |1.9| ) 
to the corresponding ^-exponent and which can be ex- 
tracted from the finite-size behavior of the longitudinal 
structure factom at q = q 3 (M). 



C. A next-to- nearest-neighbor coupling 



A next-to-nearest-neighbor coupling 



N 



H2 = 2 ^ Si ■ S; +2 , 



(1.13) 



H(h 3 ,a) = H{h 3 ) + aH 2l 



(1.14) 



does not change the position of the soft modes q± — tt and 
q 3 (M) but, changes the associated r?i(M, a), r] 3 (M 7 a) 
exponents!! A singlet triplet gap ©sens in the free field 
case ( h 3 — 0) for a > a c — 0.241 ... .113 Note, however, that 
( 1.13 ) is translation invariant and therefore the ground 
state degeneracy with momenta p = 0, tt - predicted by 
the LSM construction - still holds, i.e. the singlet ground 
state is still twofold degenerate in the singlet sector. 



D. A staggered dimer field 

A plateau in the magnetization curve at M = 1/4 has 
been found in the Hamiltonian (1.4) with an additional 



next-ptcwuearest-neighbor coupling and a staggered dimer 

fieid : ye 



H(/i 3 , a, S) = H(h 3 ) + «H 2 + SD(tt). 
The dimer operator is defined as: 



N 



D(q) = 2Y e^Si ■ S 



1=1 



i+i- 



(1.15) 



(1.16) 



Such a Hamiltonian only commutes with T 2 and there- 
fore reduces the degeneracy of the ground state. At 
M = the twofold degeneracy of the ground state is 
lifted and a gap opens between the energies of the ground 
state and the excited states.1121 

At M = 1/4 a gap opens between the ground state |0) 
and one further state, which is even under T 2 . The LSM 
construction yields a second state |1) = U|0) degenerate 
with the ground state |0), which is odd under T 2 . 



E. A periodic dimer field 

A plateauJn the magnetization curve at M = 1/6 has 
been foundll3 for a Hamiltonian of the type ( |l.4|) with a 
dimer field D(<?) of period q = 27r/3: 



D(g) = [D(g)+D(-g)]/2. 



(1.17) 



The Hamiltonian used in Ref. |IfJ can be reformulated 
as a single spin- 1/2 chain with ferromagnetic nearest 
neighbor coupling being strongly disturbed by an anti- 
ferromagnetic next-to-nearest neighbor coupling with a 
period of q = 27r/3. 

Note, that the periodicity q = 2n/3 coincides with the 
soft mode momentum q 3 (M = 1/6). Such a coincidence 
occurs in both examples pTa and IE and we conclude 



1=1 



that the special type of the periodic perturbation (1.11) 
and (1.17) is not relevant for the formation of a magne- 
tization plateau. 



2 



The situation in example 1 D is different. Here the peri- 
odicity (q = 7r) of the staggered dimer fiel d coi ncides with 
the second soft mode 2q 3 (M = 1/4) = tt [ (1.7) for k = 2], 



(A|S||jfe) 



predicted by the LSM construction. Note, however, that 
the magnetization plateau at M = 1/4 is only visible if 
the parameters a and 5 in ( 1.15 ) are appropriately cho- 
sen. In particular the magnetization plateau at M = 1/4 
seems to be absent if the next-to-nearest-neighbor cou- 
pling a is switched off. Therefore, we conclude that the 
coincidence of the periodicity q in the perturbation op- 
erator with the momentum of one LSM-soft mode is a 
necessary - but not sufficient - condition for the forma- 
tion of a plateau. 

It is the purpose of this paper to investigate in more 
detail the mechanism for the formation of gaps and 
magnetization plat eaus by means of periodic dimer- 
perturbations ( 1.17 ) of strength 2S q . The ^-evolution 
of the energy eigenvalues and transition matrix elements 
of the perturbation operator is given by a closed set of 
differential equations, which we have discussed in Refs. 
[ToUTl} The initial values for these evolution equations 
are given by the energy eigenvalues and transition am- 
plitudes for the unperturbed case (8 q = 0). 

The outline of the paper is as follows. In Sec. [n] we 
complete the discussion of the quantum numbers of the 
LSM state |fc) by investigating their S T content, where 
St is the total spin operator. We are in particular inter- 
ested in the question, whether or not the LSM state |1) 
at M = with momentum pi = po + n contains a triplet 
(S =1) or higher spin component [S T = S(S + 1)]. Sec- 
tion |n is devoted to an analysis of the LSM soft modes 
in the dimer-dimer structure factor. This analysis is used 
to fix the above mentioned initial conditions for the evo- 
lution equations. In Sec. IV we then present numerical 



results on the formation of gaps and plateaus by means 
of the periodic dimer perturbations. 

The occurrence of magnetization plateaus in spin lad- 
ders is discussed in Sec. |V|. 



II. THE LIEB, SCHULTZ, MATTIS (LSM) 
CONSTRUCTION AND THE QUANTUM 
NUMBERS OF THE DEGENERATE GROUND 
STATES. 

It has been pointed out in the introduction that the 
quantum number analysis of the states \k) = U fc |0) in 
the LSM construction is crucial to decide whether these 
states are new, i.e. orthogonal to the gro und state |0), 
or not. The transformation beha vior Jl.q ) under trans- 
lations T yields the mom enta pk (1-7) of these states. 

The operator U (1.1) obviously commutes with the 



total spin in 3-direction Sj.- Therefore all the states 
|fc) = U fe |0) have the sam e tota l Spin S T in 3-direction. 
The Hamiltonian of type (1.14) with isotropic couplings 



commutes with S T . One might ask for the S T content of 
the states \k). To answer this question we compute the 
expectation value 



(0|S X |0> = 
2N[Si(q 



(0|U tfc S T U fe 



|0) - <0|S£|0) 
k2ir/N,M) - 5i(0,M)], (2.1) 



using the considerations developed by LSM to show the 
vanishing of the energy difference (1.3). The right-hand 
side of ( |2.lD is determined by the transverse structure 
factor, exposed to an external field h 3 with magnetization 
M(h 3 ): 



Si(q,M) 



N 

E 



(2.2) 



which has been studied on fi nite systems in Ref. [It] for 
the nearest neighbor model (1.4) and in Ref. [l8] for the 
model with next-to- nearest-neighbor couplings a. 

From these investigations we conclude for the thermo- 
dynamical limit of the difference appearing on the right- 
hand side of (2.1): 



2N[Si (k2n/N, M) - Si (0, M)\ 



N- 



k 



(2.3) 



The exponent ft — ft [M, a) turns out to be zero for 
M = and a = [cf. Fig. 5(d) in Ref. [P?]]. This means, 
that the right-hand side of (^J) is non vanishing. For 
M = 0, a = the ground state is a singlet state [S T = 
S(S + 1) = 0] and (2.1) tells us that the soft mode state 
\k = 1) = U|0) with momentum it contains triplet [S T = 
S(S + 1) = 2] and higher spin components, i.e. the LSM 
construction together with (2.2) and j3k{M = 0, a = 0) = 
forbids a singlet triplet gap. 

The exponent (3k(M, a) is larger than zero for M > 
[cf. Inse t of Fig. 5(b) in Ref. [t]] In this case the right-hand 
side of (2.1) vanishes and the soft mode states \k) = 
U fc |0) have the same total spin S T = S(S + 1) as the 
ground state. 

Switching on the next-to-nearest-neighbor coupling a 
the free field exponent must be larger than zero 



/3i (0, a) > for a > a c = 0.241 



(2.4) 



since the diper phase a > a c is characterized by a singlet 
triplet gap.E3 In other words, for a > a c the degenerate 
LSM state |1) = U|0) with momentum p s + ir must be a 
pure singlet state as well. The exponent /?i(Q_L/2) can 
easily be calculated at the Majumdar-Ghosht2lE2l point 
a = 1/2 . Here we find ft (0,1/2) = 1. 

In Fig. |l| we have plotted the energy difference = 
Ei — Eo of the two singlet states |0) and |1) 



3 



0.06 



0.03 



-0.03 



-0.06 



8 


10 








12 












14 



0.5 



0.6 



0.7 0.8 
a 



0.9 



FIG. 1. The energy differences (2.5) for system sizes 
N = 8, 10, . . . , 16 and S = 0. 



As was pointed out in the introduction, a fourfold de- 
generacy of the states |fe) = U fc |0), (fc = 0,1,2,3) is 
predicted at M = 1/4. All these states have the same to- 
tal spin squared (S = N/4) and momenta p = 0, ±7r/2, tt. 
On finite systems one again observes oscillations in the 
energy differences 

lu p (S, a) = E(p, S = N/4, a) - E(0, S = N/4, a), (2.8) 

for p = 7r/2, it, if we switch on the next-to-nearest- 
neighbor coupling a as is shown in Fig. |^. For a values 
large enough these oscillations die out and the ground 
state momentum p s ,S = N/4 is supposed to be 



TV = 8, 24, 40,. 
N = 12,28,44, 
N = 16,32,48, 
N = 20,36,52, 



(2.9) 



u v (S = 0, a) = E(p s + tt, S = 0, a) - E(p s , S = 0, a), 

(2.5) 

versus the coupling a. This is an oscillating function for 
a > 0.5 with zeroes at: 

a = ai(N) < a 2 (N) < ... < a z (N), (2.6) 

where numerical data suggest that the total number of 
zeroes is given by 



llTV : N = 8, 12, 16, . 
4 | AT -2 : N = 10, 14, 18, 



(2.7) 



In the thermodynamical limit we have a dense distribu- 
tion of zeroes. The height of the maxima and minima in 
between converges to zero - a signal for the degeneracy of 
the two states |0) and |1) = U|0) in the thermodynamical 
limit as predicted by the LSM construction. 



III. SOFT MODES IN THE DIMER DIMER 
STRUCTURE FACTOR 

According to the LSM construction for the transla- 
tion invariant models with nearest and next-to-nearest- 



neighbor couplings (1.14), H(/i3,a), we expect the dis- 
persion curve 

Wq(S,a) = E(p s + q,S,a) - E(p s ,S,a), (3.1) 



to develop zeroes 



uM(/u 1 Q)=u,(S ) a) 1 



(3.2) 



at the soft mode momenta q = q^ k \M) = fe7r(l — 2M). 
If in the thermodynamical limit the scaled energy differ- 
ences: 

CP\M,a)= lim Nuj^(h 3 ,a). (3.3) 

N— >oo 



and 



s 



0.1 



\ (16,71) 
(1 6,71/2) '"\ 




FIG. 2. The energy differences (2.8) for system sizes 
N = 8, 12, 16 and momentum p = n, n/2 and S = N/4. Same 
line types denote the same system sizes. 



v{M,a) = lim N[E(p s + 2ir/N, S,a) - E(p s , S,a)}, 

TV— »oo 



are finite and non vanishing, the ratios 

ttv(M, a) 



(3.4) 



(3.5) 



yield the exponents r)( k '(M, a), which govern the critical 
behavior of the dimer dimer structure factor: 

S DD (q,M) = ^(p s ,S\B c (q)-Dl( q )\ Ps ,S). (3.6) 



Here 



D c (?) = D(?) - <J,„D(0) 



4 



is the connected part of the dimer operator ( 1.1 6| ) . For 
AI = a zero should occur at q = it and for M = 1/4 
we should find two zeroes at q — tt/2 and q — tt. This is 
indeed the case as can be seen from Fig. || where we have 
plotted the dispersion curves for M = 1/4, a = 0, 1/4. 




FIG. 3. The dispersion curve of the Hamiltonian (1.14) for 
finite system (N = 24, 28) and magnetization M — 1/4. 

Note, that the dimer operator commutes with the to- 
tal Spin and the dispersion curve (3d) describes the 
lowest lying excitations contributing to SoD{q, M, a). 

T he o -dependence of the dimer-dimer structure fac- 
tor (3.6) at M = 1/4 and a = 0,1/4 is shown in 
Fig. I|r^ A pronounced peak is found at the first soft 
mode q — tt(1 — 2M) = tt/2. The finite-size behavior 
of the structure factor at the first soft mode 



S DD (q 3 (M),M,a) ~ N i-n 3 (M,<*) 



(3.7) 



is shown for a = 0, 1/4 and M — 1/4 in the inset of 
Fig. 0. It is well described by an exponent 



7 7 W(l/4,a)=r 73 (l/4,a) 




(3.8) 



identical with the exponent 773 (M , a) in the longitudinal 
structure factor. The latter has been calculated exactly 
in the model with nearest neighbor couplings (a = 0) by 
means of Bethe ansatz solutions for the energ y eig enval- 
ues, which enter in the differences (3.1) and (3.4). The 
resulting curve 773 (M, a = 0) is shown in Fig. 2 of Ref. [|. 
The a-dependence of 773(1/4,0;) has been calculated on 
small systems (N < 28) and can be seen in Fig. 4(b) of 
Ref. g. 

According to its definition ( |3.6| ), the dimer dimer struc- 
ture factor can be represented: 



>dd \ 



( ? ,M,a) = i£Kn|D c ( g )|0>| s 



(3.9) 



in terms of transition amplitudes from the ground state 
|0) = \ps,S) to excited states |n) with momenta p n — 



p s + q and total spin S. The peak in Fig.[| tells us that 
the transition matrix elements at the first soft mode q = 
tt(1 - 2M) 



(n|D(7r(l-2Af))|0) 



JV- 



(3.10) 



with K3 = 1 — ?73(M(/i3), a)/2, diverge with the system 
size. 
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FIG. 4. The dimer dimer structure factor (B.£) at M = 1/4 
for = 12, 16, . . . , 28. The insets show the iV-dependence of 
the structure factor Sdd at q — tt/2, versus iV 1- '' 3 , with 
exponent 773 given by Eq. ( J3.8| ) . For a = 1/4 the structure 
factor diverges and for a = it is finite. The extrapolated 
value (N — > 00) is marked with a solid symbol (■). 

The fact, that there is no peak at the second soft mode, 
indicates that the corresponding transition matrix ele- 
ments (n\T)(q — 2ir(l — 2M))|0) are small - at least for 
next-to-nearest-neighbor couplings a < 1/4. The mag- 
nitude of these transition matrix elements is crucial for 
the formation of gaps and magnetization plateaus with a 
periodic perturbation T)(q). 



IV. PERIODIC DIMER PERTURBATIONS AND 
THE FORMATION OF GAPS AND PLATEAUS 

In this section we will study the impact of a periodic 
dimer perturbation D(q) (1.17): 



3, a, 



H(h 3 ,a)+26 q T>(q), 



(4.1) 



on the soft modes. We will restrict ourselves to the model 
with nearest neighbor and next-to- nearest-neighbor cou- 
plings and follow the procedure of Refs. fojll There it 



was shown that the 5 g -evolution of the energy eigenvalues 
and transition matrix elements (n|D(g)|0) is described 
by a closed set of differential equations which possess fi- 
nite size scaling solutions in the limit N — > 00, 5 q — > 0, 
x = S q N fixed. 
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This statement holds, if the periodicity q of the dimer 
perturbation coincides with a soft mode momentum: 



q = g (fe) (A/) = Ml - 2M )- 



does not vanish in the thermodynamical limit, i.e. if 
the ground state energies E(p s i,S',a,6 q ) for S 1 = S — 
1, S, S + 1 evolve in a different manner under the pertur- 
(4-2) bation<5 g . This happens exactly if the periodicity q of the 



Then a gap in the energy difference ( |3.l| ) between the 
ground state and lowest state which can be reached with 
the perturbation t)(q) at q^ k '{M) is predicted: 

w<*> (h 3 , a,S q ) = Sf fi< fc > (M, a, x). (4.3) 

The gap opens with an exponent = e^ k \h3,a), re- 
lated to the corresponding ?/-exponent (pTq) via 



(h 3 ,a) =2{4- V ^(M(h 3 ),a)\ 



(4.4) 



A test of the prediction (O) is given in Fig. 0, where we 
plotted the gap ratio 



u^(h 3 ,a,S q ) 
u>M{h 3 ,a,0) 



l + e^(h 3 ,a,x), (4.5) 



for M — 1/4, i.e. q = tt/2 and a = 0,1/4 versus the 
scaling variable x 2 ^ with the exponents 



(1/4, a) 



0.81(1) :a = 
0.64(3) :a = l/4 



(4.6) 



to 
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FIG. 5. The gap ratio (4.5) v ersu s the scaling variable 
(7v<5 E ) 2/e_1 , with e given by Eq. Q). The solid lines rep- 
resent linear fits to the small a;-behavior. 



Note that the gap ratio (4^5) is linear in x 2 ^ for small 



values of predicteeLby the evolution equations for 

the scaling functions Ej'Eil 

Let us nex t study the influence of the periodic pertur- 
bation ( 1.17 ) on the magnetization curve M = M(h 3 ). 
A plateau in the magnetization curve with an upper and 
lower critical field h 3 and h\ will emerge if 

A(5 g ,h 3 ) = h 3 l -h 3 = lim [Efr^S+l^SJ 

N^oo 

+E{p s _ 1 ,S - 1, a, S q ) - 2E{ Ps , S, a, S q )} , (4.7) 



dimer pe rturbation coincides with a soft mode momen- 
tum (fL2\). Then the ground state energy at the 'critical' 
magnetization M — Sj^/N is lowered stronger than at 
the neighboring magnetizations M = (S^ ± 1)/N. 

In Fig. ^ we show the <5 q -evolution of the magnetization 
curves for q = tt/2 and a = 1/4. The emergence of the 
predicted plateau at M = 1/4 is clearly visible. 




FIG. 6. The magnetization curve of H(h:j,a, 8 q ) for 
a — 0.25 and different <5 ? -values, determined from finite sys- 
tem sizes N = 8, 12, 16, 20. 



The scal ing _ behavior S q with e 



eW(/i3,a) of the 

difference (4/7) is governed by the critical exponent 
r}^ 1 \h 3 ,a) at h 3 = h 3 (M = 1/4) of the unperturbed 
model at the first soft mode q 3 (M) as can be seen in 
Fig. 0. 

We also looked for the (5-evolution of the magnetization 
curves for a dimer perturbation with periodicity q = ir. 
There is a plateau for M = - corresponding to the gap 
above the ground state discussed in Ref. |^. However, no 
plateau is visible at M — 1/4 for small perturbations 
5D(tt). This corresponds to the observation that the sec- 
ond soft mode at M = 1/4, q = w, does not produce any 
signature in the dimer dimer structure factor. 

These statements only hold for small perturbations 
5T)(tv). It is indeed known from Refs. |lj,|l^ that a plateau 



6 



at M = 1/4 can be enforced with a large perturbation of 
order 1. 




FIG. 7. The evolution of the difference (|OJ between the 
upper and lower critical field at the plateau M — 1/4. The 
dotted and dashed line show a fit, proportional to 8 e q with 
e = e' 1 ) given by Eq. (4.6), to small values of the external 
field 5„. 




FIG. 8. The magnetization curve of Hamiltonian ( 1. 1 ) with 
a perturbation D(27r/3) + D(-7r/3) for a — and different 
J 3 -values, determined from finite system sizes N = 6, 12, 18. 



q = 2ir/3 and q — ir/3. We found clear evidence for 
the expected magnetization plateaus at M —1/6 and 
M = 1/3. The magnetization curve for a Hamiltonian 
with a superposition of both perturbations 



D(2tt/3) + D(tt/3), 



(4.8) 



is shown in Fig. ^. Here, we find two plateaus in the 
magnetization curve at M — 1/6 and M = 1/3. 



V. MAGNETIZATION PLATEAUS IN SPIN 
LADDERS 

All the considerations we made so far for the ID Hamil- 



tonian (4.1) with nearest and next-to-nearest-neighbor 
couplings can be extended to the case 



B.(h 3 ,a l ) = B.(h 3 ) + a l li l , 



(5.1) 



where we substitute the next-to-nearest-neighbor cou- 
pling by couplings over / lattice spacings 



N 



(5.2) 



For I finite the position q 3 (M) of the first soft mode - 
generated by the LSM construction - will not change, 
the corresponding ^-exponent might change. According 
to our experience with the case I = 2, we expect that a 
slightly frustrating coupling enhances the singularity in 
the dimer dimer st ruct ure factor at q — q 3 {M). Hamilto- 
nians of the type (5.1) are interesting, since they can be 
viewed as a spin ladder system with I legs, as is shown in 
Fig. 0. 
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FIG . 9 . A spin ladder system described by Hamilto- 
nian (5.1) with / legs, and additional diagonal couplings 
(dashed lines). The coupling strength between the legs is 
one unit. 



We have c omp uted magnetization curves for the 
Hamiltonian (|4.1|) with perturbations D(g) of period 
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They differ, however, from usual spin ladder systems 
with I legs, owing to the diagonal (dashed) couplings [(I <-> 
I + 1), (2/<->2i + 1), . . .], which are inferred by the helical 
boundary conditions. 

Indeed, these additional couplings change the physical 
properties. Spin ladder systems with helical boundary 
conditions are gapless - irrespective of the number of 
legs. This statement holds if the couplings a/ along the 
legs is chosen properly, e.g. the two leg system with I = 2 
is gapless for a < a c = 0.241 .... Conventional ladder 
systems are known to be gapless for an odd number of 
legs, but to be gapped for an even number of legs. 

This fundamental difference becomes clear, when we 



add a special dimer field to (5.1) 



jv 



n+l, 



(5.3) 



71=1 



which only affects the diagonal couplings [(/ *-* I + 
l),(2J~2?+l),...],inFig.g 



J n+l ~ J n ■ 




.1-1 



(5.4) 
(5.5) 



The periodicity (5.5) of the couplings J^p is given by a 
Fourier series 

[1/2] 

J« = cos(2Tmj/l)5(q = 2nj/l), (5.6) 

3=0 

where [1/2] is the largest integer sm aller than 1/2. The 
Fourier coefficients 5(q) follow from (5.4), e.g. for I = 2 
we find: 



S(q = 0) = 5{q = tt) = 8/2. 
The dimer perturbation: 

<5D( 2 ) = ^[D(0)+DW] 



(5.7) 



(5.8) 



generates a plateau at M = (1 — q/n)/2 = 0. Therefore, 
the gap - typical for the two leg ladder - appears imme- 
diately if we switch on the dimer field (5.3), which breaks 



the translation invariance of the ID system. 

Let us next consider the three leg ladder (1 = 3). The 
Fourier coefficients turn out to be 



5(q = 0) 



C3 



-6, 6{q = 2tt/3) 



1 



1 + C 3 



S, (5.9) 



with C3 = cos( 7r/3). The q = 2n/3 component in the 



1 + C 3 
COS (V/ 

dimer field (plj|): 

5T>( 3) = 5{q = 0)D(0) + 5{q = 27r/3)D(27r/3), (5.10) 



generates a plateau at M = 1/6. This is exactly the 
plateau found in Ref. Efi. 



Note, that the Fourier decomposition of SD^ contains 
in general a q = n component for even I and no q = tt 
component for odd I. We therefore find a gap at M = 
for ladders with an even number of legs but no gap for 
ladders with an odd number of legs. 

In summary we c an say, the Fourier decomposition of 
the dimer field (O) 



<5D« 



[1/2] 

3=0 



2nj/l)B(q = 2-Kj/l), 



(5.11) 



tells us, where to expect plateaus in the magnetization 
curve of a spin ladder system with I legs. The position 
of plateaus can be seen in Table |. 

The number of plateaus increases with the number 
of legs and one is tempted to suggest that in the two 
dimensional limit I — > 00, the magnetization curve is 
again a continuous function. It should be noted, how- 
ever, that the LSM-construction with the operator (|l . l| ) 
breaks down in the combined limit N — * 00, k = vTV in 



the sense that (1.3) does not hold. 

The magnetic properties of the two dimensional 
Heisenberg model with helical boundary conditions at 
T = have been studied in Ref. ^2|. Concerning the 
isotropic model with nearest neighbour coupling, there is 
no indication for a plateau in the magnetization curve. 

Finally, let us mention that there is a second way to 
map ladder systems with I legs onto one-dimensional sys- 
tems with far reaching couplings: 

H(/ 1 3,t ; ) = H(/ 13 ) + t ; H 7V/; . (5.12) 

The couplings for the I leg system are shown in Fig. [lO], 

12 3 N/l-1 N/l 
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1 1 
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N-I+l N-l+2 N-l N 

FIG. 10. Spin ladde r syst em with / legs described by Hamil- 

tonian H(h 3 ,Ti) [Eq. (|5.12|)]. 



and for the two leg system in Fig. pj| , 

12 3 N/2-1 N/2 



N/2+1 N/2+2 N 
FIG. 11. Spin ladder system with two legs described by 
Hamiltonian H(/i3,T2). 
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The latter differs from the conventional two leg sys- 
tem (with periodic boundary-conditions) only by a twist 
at the boundary, which should not change the physical 
properties in the thermodynamical limit. Therefore, we 
expect a gap in this system. The appearance of a gap 
in the systems with an even number I of le gs or iginates 
from the second term in the Hamiltonian (|5 . 1 2 ) . If we 



repeat the calcul ation of the expectation values (1.3) for 
the Hamiltonian ( 5.1 2| ) we find: 



(fc|H|fc) - (0|H|0) = OiN- 1 ) 

N 

+ ^((o|u fc s„s„ +jV/; ut' i _s„s 



nJ n +N/l\ 



N 



71=1 



The coefficient /* 
unless 



[cos(27rfc/Z) — 1]tj does not vanish 



k = 1,21,.. 



(5.14) 



This means in particular that the LSM-operators U fe , 
k = 1, ... ,1—1 do not create states, which are degenerate 
with the ground state in the thermodynamical limit. 

Those situations, where the ground state degeneracy is 
lifted completely, are of special interest. They occur if the 
momenta of the states \l) = U z |0) [pi = ln(l - 2M) + p s ] 
and of the ground state |0) (p s ) differ by a multiple of 
2-7T, i.e. for 



1(1 

2 K 



2M) e Z. 



(5.15) 



The condition ( 5.15 ) is satisfied exactly for the I and M 
values, listed in Table |, i.e. for those values, where we 
expect a plateau in the magnetization curve. 



VI. CONCLUSION AND PERSPECTIVES 

In this paper we tried to elucidate the mechanism 
which generates gaps and plateaus in spin 1/2 antifcr- 
romagnetic Heisenberg models with nearest and next 
to nearest neighbor couplings of strength a. A priori 
these models have no gap in the presence of a homo- 
geneous field /13 > h^(a) above the critical field h^(a), 
which is needed to surmount the singlet triplet gap for 
a > a c — 0.241 .... The LSM construction predicts the 
existence of soft modes (zero energy excitations) at wave 
vectors qW(M) = Wl - 2M), k = 1, 2, 3, . . . . 

It is shown in Sec. [jj that the total spin squared = 
S(S + 1) of the lowest excited states at the soft mode 
momenta is the same as that of the ground state for S — 
MN. The soft modes are therefore expected to generate 
signatures in the dimer dimer structure factor, since the 
dimer operator does not change the total spin squared. 



Indeed, for M = 1/4 a pronounced peak is seen at the 
first soft mode q = gW(l/4) = tt/2 if a < 1/4, indicating 
a large transition matrix element (l|D(7r/2)|0) between 
the ground state - with momentum p s - and the first 
excited state with momentum p s + n/2. 

The second soft mode q = g( 2 )(l/4) = tt, however, 
does not produce any visible structure in the dimer dimer 
structure factor (for a < 0.25). Here the relevant tran- 
sition matrix elements (l|D(7r)|0) between the ground 
state |0) and the first excited state with momentum p s +n 
are small. The situation is different for large next-to- 
nearest-neighbor couplings a. 

The magnitude of the transition matrix elements 
(l|D(q)|0) is crucial for the efficiency of the mechanism 
to generate a plateau in the magnetization curve at a 
rational value of M by means of a periodic perturba- 
tion 5qD{q). According to the criterium of Oshikawa, 
Yamanaka and AffleckO a plateau at M is possible if q 
meets one of the soft modes q = q^ k \M), k = 1, 2,3, 

Our numerical analysis shows that the width of the 
plateau - i.e. the difference of the upper and lower criti- 
cal field (4.7) - depends on the magnitude of the transi- 
tion matrix element in the unperturbed model (S q = 0). 
Indeed these matrix elements enter as initial conditions 
in the differential equations [(2.4), (2. 5)] and [(2. 2), (2. 3)] 
in Refs. [l(||ll], which describe the evolution of gaps and 
plateaus under the influence of a periodic perturbation 
5 q T)(q). We have also demonstrated that a superposition 
of two periodic perturbations [D(27r/3) + D(-7r/3)] gener- 
ates two plateaus in the magnetization curve exactly at 
those magnetization values (M = 1/6, 1/3) where the pe- 
riod (q = 27r/3,7r/3) in the perturbation coincides with 
the first soft mode momentum q — q^(Al). 

Ladder systems with I legs [cf. Fig. can be inter- 
preted as one-dimensional systems with additional cou- 
pling s ov er I l attic e spacings and a dimerized perturba- 
tion (5.3) and (5.4). The latter breaks translatio n inv ari- 
ance of the ID system and t he F ourier analysis (5.11) of 
the dimerized perturbation (5^) reveals the occurrence 
of magnetization plateaus [cf. Tab. |] in spin ladder sys- 
tems with I legs. 

In this paper we concentrated our investigations on 
different spin- 1/2 systems that all had in common that 
they reduce in the unperturbed case to critical Heisen- 
berg chains. It should be pointed out that there exist 
exact solutions on other (multichain spin- 1/2 antiferro- 
magnetic Heisenberg-like) modelsj23 where the existence 
of magnetization plateaus is still unclear. We remark 
that the application of the method we presented here, is 
not limited to the cases we discussed. However, besides 
looking for the existence of soft modes it turns out to be 
needed to discuss the strength of the transition matrix 
elements. Here, we cannot offer general prescription and 
therefore the dynamics of each different model of interest 
has to be treated separately. 
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